Multipole moments of bumpy black holes 
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General relativity predicts the existence of black holes, compact objects whose spacetimes depend 
on only their mass, spin, and charge in vacuum (the "no hair" theorem). As various observations 
probe deeper into the strong fields of black hole candidates, it is becoming possible to test this 
prediction. Previous work suggested that such tests can be performed by measuring whether the 
multipolar structure of black hole candidates has the form that general relativity demands, and 
introduced a family of "bumpy black hole" spacetimes to be used for making these measurements. 
These spacetimes have generalized multipoles, where the deviation from the Kerr metric depends on 
the spacetime's "bumpiness." In this paper, we show how to compute the Geroch-Hansen moments 
of a bumpy black hole, demonstrating that there is a clean mapping between the deviations used 
in the bumpy black hole formalism and the Geroch-Hansen moments. We also extend our previous 
results to define bumpy black holes whose current moments, analogous to magnetic moments of 
electrodynamics, deviate from the canonical Kerr value. 

PACS numbers: 04.25.Nx, 04.70.Bw 



I. INTRODUCTION 



In many areas of physics, multipolar expansions are used as tools for describing the shape of a distribution of 
<sj i matter or energy, or for describing the behavior of a potential function. Multipole moments are most commonly used 
^ ' to describe fields whose governing equations are linear, since the functions describing the angular behavior are typically 
eigenfunctions of the angular piece of the governing differential operator. For example, the Newtonian gravitational 
potential $ arising from a matter distribution p must satisfy Poisson's equation: 



AirGp (interior), 



>: v'* = | ; (i.i) 

| In the exterior region, we can write $ as a sum over multipolar contributions: 
00- *(r,M) = -G^ jj+i ■ (1-2) 

o . 

The coefficients Mi m are mass multipole moments. By matching the expansion of $ on the boundary to a similar 
expansion for the interior, they can be shown to describe the angular distribution of the mass of the source. (Through- 
out this paper, we will restrict ourselves to axisymmetric spacetimes, for which the axial index m must be zero; we 
ignore it in what follows.) 

Using multipole moments to describe gravity in general relativity is not as simple, thanks largely to the nonlinear 
nature of the governing equations. Geroch [l[ and Hansen [2j developed a very useful multipolar description for 
spacetimes of isolated, stationary, axisymmetric objects in GR in terms of scalar multipoles. Their definition applies to 
spacetimes that are asymptotically flat; for such spacetimes there is a well-defined "large r" limit in which multipoles 
can be defined in a way that roughly accords with our usual intuition. When one computes the Geroch-Hansen 
moments of a source, one finds that its spacetime is described by a family of mass moments Mi , very similar to those 
appearing in our Newtonian expansions in Eq. (|1.2j) . as well as a family of current moments Si. For a fluid body, the 
current moments describe how the matter flow is distributed through the compact body, much as magnetic moments 
describe how electric current is distributed through an electromagnetic source. These moments can be conveniently 
combined in the complex moment Aii = Mi + iSi. 

For a generic source, the moments M.i are unconstrained. 1 As we will discuss in more detail in the following section, 
for a Kerr black hole the moments take a particularly simple form: 

Mi = M(ia) 1 , (1.3) 



If we assume the spacetime is reflection symmetric, then only the even mass moments and odd spin moments can be nonzero. In this 
paper, we do not restrict ourselves to reflection symmetric spacetimes. 
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where M.\ is the Zth moment, M is the total mass of the black hole, a is its spin parameter, and we are using units 
where G = c = 1. (We neglect the astrophysically uninteresting possibility of a black hole with macroscopic charge.) 
This is a statement of the "no-hair" theorem: the spacetime of a Kerr black hole is completely described by its mass 
and spin Q. 

The sharply constrained nature of the Kerr multipoles implied by Eq. (( 1 . 3|) suggests that this relation may be 
useful as a test of black hole spacetimes: if the spacetime is Kerr, then knowledge of only two moments is needed to 
determine all of the others. Ryan [J] was the first to build a scheme to test this idea, constructing spacetimes in which 
all of the moments were arbitrary. This scheme was sufficient to prove the principle of the idea, but did not work 
well for building spacetimes good deep into the strong field. Collins and Hughes [5| noted that, if general relativity 
correctly describes black hole candidates, then testing their nature amounts to trying to falsify the hypothesis that 
they are Kerr black holes. They suggested formulating black hole tests as a null experiment by examining spacetimes 
for which 

Mi = M(ia) 1 + SMi , (1.4) 

and using measurements to test whether 6A4i — 0, as it should if they are Kerr black holes. 

Collins and Hughes formulated a simple version of this test by showing how to deform the multipoles of a non- 
rotating black hole. Glampedakis and Babak Q extended this idea by modifying the quadrupole moment of a Kerr 
black hole. Though not as generic in the form of the moment that can be modified, this significantly improved 
the astrophysical relevance of this test by allowing for spacetimes with angular momentum. Vigeland and Hughes 
[?J then improved the Collins and Hughes construction, using the Newman- Janis algorithm to transform bumpy 
Schwarzschild black holes into bumpy Kerr black holes. Their approach allows one to vary any of the spacetime's 
mass moments, an important consideration since black hole candidates may agree with the Kerr metric's multipoles 
up to some ? max but differ for I > l max . 

We have chosen to focus on the Geroch-Hansen moments, as opposed to the tensor moments described by Thorne [§], 
for the sake of simplicity. However, as we discuss in Sec. HH Geroch-Hansen moments are only defined for spacetimes 
for which R^ v = 0. As we show in Sec. IIVI and I VI) when we add bumps to the Kerr metric, the spacetime is no 
longer vacuum in GR. It will turn out that for the perturbations we have considered, the nonvacuum nature is of 
sufficiently "high order" (in a sense made precise in Sec. IIV|) that we can still define Geroch-Hansen moments up to 
some order. In order to study the multipolar structure of other kinds of perturbations to the Kerr metric, we would 
need to calculate the spacetime's tensor moments. For example, this would facilitate analyzing the "bumpiness" of 
black holes in other theories of gravity, such as dynamical Chern-Simons extensions to general relativity. We leave 
this for future work. 

The purpose of this paper is twofold. First, we map the mass perturbations described in [?J to perturbations in the 
Geroch-Hansen moments SAii. In constructing a bumpy black hole, we work in a few different coordinate systems, 
and specify the spacetime's bumps by adding a function which is proportional to a spherical harmonic in one of 
those coordinate systems. Our goal is to demonstrate that this choice maps in a natural way to the spacetime's 
Geroch-Hansen moments. For example, a bump which is built from an I — 2 spherical harmonic changes the even 
Geroch-Hansen moments above the / = 2 moment. In principle, one could construct a change to a single Geroch- 
Hansen mass moment by including multiple appropriately weighted spherical harmonic terms in the bumpy black hole 
spacetime. 

Second, we show how one can build bumpy black hole spacetimes for which the current moments are modified. 
This expands the domain of tests that one can perform, making it possible to test the full range of the gravitational 
moments. For example, Yunes and Pretorius [Toj showed that dynamical Chern-Simons extensions to general relativity 
lead to a rotating black hole solution that is not the Kerr metric; in the slow rotation limit, the spacetime's moments 
are the same as Kerr for I < 3, but its I — 4 current moment differs from the Kerr value (Tlj . Having control over 
both a spacetime's mass and spin moments completes our ability to build bumpy black holes and allows us to test 
the Kerr metric in great detail. We show that the spin moment perturbations behave similarly to the mass moment 
perturbations; if we add an order I spin moment perturbation, which we define in Sec. [V] and I VII we leave the spin 
moments unchanged up to order I. When we allow for both mass moment and spin moment perturbations, the bumpy 
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Kerr spacetime is described by three perturbation potentials ip\, 71 , and a\ : 



ds 2 



(1 + 2Vi) 1 



2Mr 



dt 2 - lx - 



Aa 2 Mr sin 2 6 



AS 



A 



AaMr sin 2 6> / r 2 +a 2 2Mr A 
1 + 2-0i - 71) h 4(7i J- 1 



eft 



+(l + 2 7 i -2V-i)^-(ir 2 +2 7l 



2Mr(r 2 + a 2 ) 



(r 2 + a 2 ) 2 - a 2 A sin 2 6 + (71 - tpi) 



AE 



ga 2 M 2 r 2 gin 2 q 



a sin 2 6 dr d<j> + (1 + 2 7i - 2^i)S d6> 2 



AaMr ( . 2Mr(r 2 + a 2 ) 
-di A+ ^— '- 



£ - 2Mr 



sin 2 6» 



£ - 2Mr 
2 



-201 



£ 2 A 



£ - 2Mr 



(1.5) 



where £ = r 2 + a 2 cos 2 (9, A = r 2 — 2Mr + a 2 , and we are using Boyer-Lindquist coordinates. In the following sections, 
we will treat mass perturbations and spin perturbations separately; we include the complete bumpy Kerr metric here 
for completeness. 

The remainder of this paper is organized as follows. In SecUH we review in detail how Geroch-Hansen moments are 
calculated, and demonstrate the procedure on the Kerr spacetime. We then apply this procedure to the Schwarzschild 
spacetime with mass perturbations in Sec. Mil followed by a Kerr spacetime with mass perturbations in Sec. IIVI 
Finally, we present our procedure for perturbing the spin moments of both Schwarzschild (Sec. |V]) and Kerr (Sec. IVI|) 
black holes. We conclude in Sec. I VIII by discussing some directions for future work in this problem. We adopt the 
following conventions in this paper. We work in geometrized units where G = c = 1. When writing tensors, we use a 
Greek letter to indicate a spacetime index, and a Latin letter to indicate a spatial index. 



II. COMPUTING GEROCH-HANSEN MOMENTS 

We begin our analysis by showing how to compute the Geroch-Hansen moments for the spacetime of a compact 
object. As an important example, we demonstrate the procedure on the Kerr metric, which in Boyer-Linquist 
coordinates is written 



ds 2 = -11- 



2Mr 



dt 2 



AaMr shr 



dt d(f> + —dr 2 + T,d6 2 + [(r 2 + a 2 ) 2 - a 2 A sin 2 



sin 



dcf> 2 . (2.1) 



In order to compute the Geroch-Hansen moments, the spacetime must have a timelike Killing vector and be 
asymptotically flat. Let the spacetime's Killing tensor be K a , and let the manifold V be the 3-surface orthogonal to 
this vector. The metric on V can be written 



Xgij + KiKj 



(2.2) 



where A = —K a K a is the norm of K a . To calculate the moments, we perform a conformal transformation to map 
infinity onto a point A. The space is asymptotically flat if it can be conformally mapped to a 3-space V which satisfies 

a 

(i) V = V U A, where A is a single point; 

(ii) hij = Q 2 hij is the conformal metric; 



(iii) n| A = 0, DiQ, 



= 0, DiDjSl 



= 2h 



where f2 is the conformal factor and Di is the derivative operator associated with hij. The conformal metric has the 
form [H 



ds 2 = dr 2 + r 2 d6 2 + r 2 sin 2 9 e' 2 ^ d<b 2 



(2.3) 



The function /3(r, 9) parametrizes the deviation of the conformal metric from sphericity. 
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In order to construct the spacetime's Geroch-Hansen moments, we need its Ernst potential [HI and the conformal 
factor fl. We begin with the Ernst potential, which we build from the norm A and twist u> of K a . For the Kerr metric, 
the norm is given by 



A 



1 



2Mr 



(2.4) 



The twist is related to the "generalized curl" of the timelike Killing vector, 
From the Bianchi identities, we can write the curl as [l4[ 



V[ a W£] = —€ QL fi 1 sK 1 R V K V 



(2.5) 



(2.6) 



For a spacetime that is vacuum in GR, the condition = implies that V[ a ww = 0. This allows us to write 
w a = V q cj, where the scalar function oj is the twist of K a . (If ^ 0, we cannot construct the Ernst potential. We 
will discuss this issue in more detail in Sec. II VI and I VII ) For the Kerr spacetime, we find 



2Ma cos 9 



We then combine the norm and twist into the complex quantity 

e = A + iuj 

2Mr 2Macos9 
= 1 i 



(2.7) 

(2.8) 
(2.9) 



where on the second line we have specialized to Kerr. From this, two definitions of the Ernst potential appear in the 
literature: Ref. [13| defines it as 



1 



1 - e 



while [15] defines the Ernst potential as 



1 



1 + e 

7-1 



(2.10) 



(2.11) 



5 ) . We find the potential £ to be most useful for computing 



These definitions arc simply related to one another (£ 
multipoles of bumpy black hole spacetimes, but we use S to perturb a spacetime's current multipoles. 
Next, we must find the conformal factor J7. We begin by defining a new radial coordinate R according to 



r = R~ 



1 + MR 



M 2 - a 2 



-R- 



The point A corresponds to R = 0. The conformal factor is given by 

R 2 



M 2 -a 2 



(2.12) 



(2.13) 



R 2 ) 



a 2 R 2 sin 2 6 



and the conformal metric is given by 

ds 2 = dR 2 + R 2 d9 2 + R 2 sin 2 9 d4> 2 



1 - 



4aR sin 9 



4 - (M 2 ~ a 2 )R 2 



(2.14) 



2 Normally, the curl of a vector is only defined in a 3-dimensional vector space; we follow the lead of Ref. |15| in generalizing the notion 
of the curl here. 
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This corresponds to /3 equal to 



1 - 



AaR sin 9 



4 - (Af 2 - a 2 )i? 2 



(2.15) 



As shown by Backdahl and Herberthson [12j , the multipole moments can be computed from derivatives of a function 
y, which we now describe. Begin by defining 0, a conformally weighted variant of the Ernst potential £, 



and new cylindrical coordinates z and p: 



z = R cos 9 
p = Rsin9 . 



(2.16) 



(2.17) 
(2.18) 



We now write the potential </> as a function of these variables, 4>(z, p), and introduce yet another variation: 

<f> L {R) = 4>{R,iR) . (2.19) 
We need to define a few more functions related to the metric: 



P L (R) = p(R,iR) , 



k l (R) = -In 



1 - R 



R e 2f3 L (R') _ 1 
W 2 



dR! - RC 



Pl{R) 



(2.20) 
(2.21) 



where j3 is defined in Eq. (|2.3[) and C is the integration constant. We can choose the gauge so that C — 0. The 
multipoles are calculated from the function 



y(R) = e-«^/ 2 cj> L (R) , 



with the Zth multipole moment given by 



_ 2H\_ <fy 
1 ~ (20! dp 1 



(2.22) 



(2.23) 



p=0 



where p(R) = R e K ^'^^ R \ 

For the Kerr spacetime, the potential is given by 



MR) 



M(l + iaR) 
(1 + a 2 i? 2 ) 3/4 



The functions /3l (R) and kl (R) are given by 



(3 L {R) = -\n[l + a 2 R 2 ] 



k l (R) = -In 



2 R 2 ] 
1 + a 2 R 2 



(1 - a 2 R 2 ) 2 



Then the variable p is given by p — R{1 — a 2 R 2 ) 1 and 



y(R) = 



MVl - a 2 R 2 



1 — iaR 

In this case, the multipoles can be written compactly as 

Mi = M(ia) 1 . 



(2.24) 

(2.25) 
(2.26) 

(2.27) 
(2.28) 
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III. PERTURBATIONS TO THE MASS MOMENTS: SCHWARZSCHILD BACKGROUND 



In this section and the following one, we apply this procedure to compute the multipoles of bumpy black hole 
spacetimes as presented in Ref. 0- We begin with a bumpy Schwarzschild black hole: 



ds 2 = -(1 + 2^(1-—) 

+(1 + 271 - 2V>i) r 2 d6 2 + (1 - 2V>i) r 2 sin 2 9 d<j) 2 



2M 

dt 2 + (1 + 2 7l - 2ifa) (1 

r I \ r 



dr 2 



(3.1) 



The perturbation to the mass moments is described by the potentials ipi and 71, with rpi = 71 = corresponding to 
a Schwarzschild black hole. We further restrict ourselves to perturbations for which the metric satisfies the vacuum 
Einstein equations to first order; thus the potential ipi must satisfy Laplace's equation: 



d 2 ip! 1 dipi d 2 ipi 







dp 2 p dp dz 2 
where the Weyl coordinates (p, z) are related to Schwarzschild coordinates by 



(3.2) 



P 



rsm.9\ 1 — 



2M 



(3.3) 

z = ( r -M)cos0 . (3.4) 

Because ip\ satisfies Laplace's equation in Weyl coordinates, we can take it to be a spherical harmonic times an 
appropriate power of 1/y^p 2 + z 2 . Once we specify tpi, the potential 71 is calculated by integrating a constraint 
equation involving the perturbation ip\ and the background potentials ipo and 70, which are given by 



Mr,0) = iln(l-^ 



, 1, / M 2 sin 2 



(3.5) 
(3.6) 



See for further details. 

The norm and twist of the timelike Killing vector field, to first order in the perturbation, are given by 



The Ernst potential is given by 



A 

u = 
M 



2M 

1 + 2^1 1 

r 



The conformal factor becomes 



and the conformal metric is 



M 



n = R 2 1 - 



M 



r - M 



M 2 



-R 



(l-7i) 



ds 2 = dR 2 + R 2 d9 2 + (1 - 2 7 i) R 2 sin 2 d(j> 2 
where R has the same definition as for an unperturbed Schwarzschild spacetime, 

M 2 R 2 " 



r = R~ 



1 + MR 



(3.7) 
(3.8) 

(3.9) 

(3.10) 

(3.11) 
(3.12) 
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A. I — 2 mass perturbation 

Now we calculate the perturbations to the moments for particular choices of tpi and 71. We begin by considering 
a solution to Eq. (I3.2[) that has the form of an I = 2 spherical harmonic: 



v4 =2 (p^) 



5 2 M 3 



1 



7T (/9 2 + Z 2 ) 3 / 2 



3z 2 



pi _|_ z 2 



- 1 



We use Eq. (|3.3j) and (|3.4j) to write ^i -2 m terms of Schwarzschild coordinates: 



B 9 M 3 



3(r - M) 2 cos 2 



d(r,9y 



4 V nd{r,6) 3 
where 

d(r, 6») ee (r 2 - 2Mr + M 2 cos 2 6>) 1/2 . 
The potential 71 that corresponds to this choice of ipi is 

(r - M) [c 20 (r) + c 22 (r) cos 2 0] _ 
2 d(r,6») 5 



7i(r,0) = B 2 \j- 



where 



c 20 (r) = 2(r-M) 4 -5M 2 (r-M) 2 + 3M 4 
c 22 (r) = 5M 2 (r - M) 2 - 3M 4 . 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 
(3.18) 



We calculate the multipole moments for the bumpy Schwarzschild spacetime by following the procedure laid out in 
Sec. HU The first few multipole moments are listed below: 



5 Mo 
5 Mi 

5M 2 
SM 3 

5Mi 
SM 5 




4 



-B 2 M a 

7 V TT 



= 



(3.19) 
(3.20) 

(3.21) 
(3.22) 

(3.23) 
(3.24) 



Thus an I = 2 mass perturbation in the Weyl sector changes only the even Geroch-Hansen mass moments with I > 2, 
and the perturbations to the moments depend on magnitude of the perturbation J3 2 . 



B. / = 3 mass perturbation 



For an I = 3 perturbation in the Weyl sector, the potentials ip-y and 71 are given by 



where 



B-iM 4 7 



nd(r,e) 4 



5(r - M) 3 cos 3 9 3(r - M) cos ( 



d{r,0) 3 



d(r,6) 



b 3 m 5 Ft 



C3o{r) + c 32 (r) cos 2 9 + c 34 (r) cos 4 6 + c 36 (r) cos 
d(r,6) 7 



C30W = -3r(r-2M), 

c 32 (r) = 10r(r - 2M) + 2Af 2 

C34M = -7r(r - 2Af) , 

c 36 (r) = -2M 2 . 



(3.25) 
(3.26) 



(3.27) 
(3.28) 
(3.29) 
(3.30) 
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The first few muitipole moments are listed below: 



5Mo 


= 0, 


(3.31) 


5 Mi 


= 0, 


(3.32) 


8M 2 


= 0, 


(3.33) 


SM 3 


= -IwJl ■ 


(3.34) 


5Mi 


= 0, 


(3.35) 


SM 5 




(3.36) 



An I = 3 mass perturbation in the Weyl sector changes only the odd Geroch-Hansen mass moments with I > 3 by an 
amount proportional to B3. 



IV. PERTURBATIONS TO THE MASS MOMENTS: KERR BACKGROUND 



Now we repeat this procedure for a bumpy Kerr black hole with perturbed mass moments. We generate this 
spacetime by applying the Newman- Janis algorithm Q to the bumpy Schwarzschild spacetime, yielding [3] 



ds 2 



-(l + 2^i) 1 



2Mr\ , 9 4a 2 Mr sin 2 9 , , , N 4aMrsin 2 
- I dt 2 - 71 — dt dr - (1 + 2^i - 71) ^ 



S J 



+ (l + 2 7 i -2i)i)-dr 2 +2 7l 



1 + 



AS 

2Mr(r 2 + a 2 ) 



(r 2 + a 2 ) 2 - a A sin 2 + (71 - Vi) 



8a 2 Af 2 r 2 sin 2 
£-2Mr 



a sin 2 6> dr ^ + (1 + 2 7i - 2^i)S d6> 2 
E 2 A 



2^i 



2Mr 



sin 2 9 



(4.1) 



where, as before, A = r 2 — 2Mr + a 2 and S = r 2 + a 2 cos 2 8. This spacetime limits to the Kerr metric when ipi — > 
and 71 — ► 0; in the limit a — > 0, it reproduces the bumpy Schwarzschild metric. 

As with the bumpy Schwarzschild case, we need to compute the conformal factor and the Ernst potential for the 
spacetime. The conformal factor is 



n = r 2 



1 - 



M 2 



-R' 



- a 2 R 2 sin 2 



-1/2 



(l-7i)- 



(4.2) 



To compute the Ernst potential, we need to compute the norm and the twist of the timelike Killing vector. The norm 
is straightforward; it is given by 



A 



1-^1(1 + 2*) 



(4.3) 



Notice that the portion of the norm proportional to the perturbation falls off as r~( l+1 \ 

The twist is much trickier. Formally, it does not exist; the construction detailed in Ref. for making a bumpy 
Kerr black hole does not leave the spacetime vacuum in GR. Consider, for example, the / = 2 perturbation, described 
in more detail below. If we compute the Einstein tensor and enforce the Einstein equation = SttT^u, we find 
that, in the large r limit and to leading order in a, the spacetime has a stress-energy tensor whose only non-vanishing 
components are 



3 

8tt~ 



Tod, = -^-aB 2 M 4 \l — 



5 (3 - 5 cos 2 6) cos 9 sin * 9 
7r r 5 



(4.4) 



It may be possible to generalize Geroch-Hansen moments to nonvacuum spacetimes, but that is beyond the scope of 
this paper. For our purposes, it is sufficient to say that for the kinds of perturbations we have considered, the spacetime 
approaches vacuum very rapidly in the region where we need to use the twist. In general, since the stress-energy 
tensor is constructed from two derivatives of the metric, an order I perturbation to the metric produces non-zero 
stress-energy tensor which fall off as r~( l+3 \ This enables us to define the spacetime's twist to the order necessary to 
define the Geroch-Hansen moments to order / + 1. 
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We divide the curl into two parts: the gradient of a scalar function to' plus some correction term v a : 

u„ = V a w' + v a . (4-5) 
For an order I mass perturbation, the scalar function &/ is the same as for the unperturbed Kerr spacetime: 

2aM cos 9 



uj = — - 



(4.6) 



The correction term v a falls off as r . Since, at large r, the portion of the norm proportional to the perturbation 
falls off as r~( l+1 \ the correction to the curl can be neglected. We thus treat uj' as the spacetime's twist. 



A. 1 — 2 mass perturbation 



We generate an I = 2 mass moment perturbation on a Kerr background by applying the Newman- Janis algorithm 
to the potentials on a Schwarzschild background, Eq. (I3.14[) and (|3.16[) . This yields 



B 2 M 3 IJ 1 



= B 2 \ - 



4 V 7T d(r, 9, a) 3 
5" r 



3L(r, 9, a) 2 cos 2 9 



- 1 



d(r,9,a) 2 

L(r, 9, a) [c 2 o(V, a) + c 2 2(r, a) cos 2 9 + c 2 i{r, a) cos 4 9] 
2 d(r,6»,a) 5 



where 



and 



d(r,6,a) = Vr 2 - 2 Air + (A/ 2 + a 2 ) cos 2 9 , 
L(r,9,a) = yj(r - M) 2 + a 2 cos 2 9 , 



c 20 (r, a) = 2(r - Af) 4 - 5Af 2 (r - Af) 2 + 3A/ 4 , 

c 22 (r,a) = 5Af 2 (r- Af) 2 - 3A/ 4 + a 2 [4(r - A/) 2 - 5Af 2 ] , 

c 24 (r,a) = a 2 (2a 2 + 5Af 2 ) . 



The Ernst potential is given by 

M B 2 M 3 I5r(r-2M) 



r — M — ia cos ( 



7T (r-M) 2 



-(r - 2M) - (3r 2 - 6Mr + 2Af 2 ) cos 2 
d(r, 9) 5 



(4.7) 
(4.8) 



(4.9) 
(4.10) 



(4.11) 
(4.12) 
(4.13) 



(4.14) 



We calculate the multipole moments following the procedure outlined in the previous sections. The changes to the 
multipole moments are 



6M 
6 Mi 

SM 2 
SM 3 







(4.15) 
(4.16) 

(4.17) 
(4.18) 



The lowest order multipole that is affected is the I = 2 mass moment, and the change is the same as for an I = 2 
mass perturbation on a Schwarzschild background. The multipole moments are well-defined up to the A/1 3 moment; 
the higher order moments are not well-defined because of the presence of a fluid that behaves in the large r limit as 
an I = 4 multipole. 



10 



B. / = 3 mass perturbation 



We generate an I = 3 mass moment perturbation on a Kerr background by applying the Newman- Janis algorithm 
to the potentials on a Schwarzschild background, Eq. (I3.25|) and (|3.26[) . which yields 



^= 3 (r,l 



nd(r,9,a) 4 



5L(r, 9, a) 3 cos 3 9 3L(r,9,a)cos9 



d(r,9,a) 



B 3 M 5 Ft 



— cos( 



C3o(r, a) + c 32 (r, a) cos 2 



d(r, 9, a) 
c 34(r, a) cos 4 



C36(r, a) cos 6 9 



d(r,e,a)"' 



where 



c 30 (r,a) = -3r(r-2A0 , 

c 3 2{r,a) = 10r(r-2M) + 2M 5 -3a 2 , 

c 34 (r, a) = -7r(r - 2M) + 10a 2 , 

c 36 (r,a) = -2M 2 - 7a 2 . 



The Ernst potential is given by 
M 



B,M 4 



r — M — ia cos ( 



7 r(r - 2M) 
n r-M 



3r(r - 2M) - (5r 2 - lOMr + 2M 2 ) cos 2 



d{r,9Y' 



The changes to the multipole moments are 



6M = 

6Mi = 

SM 2 = 

6M 3 = - 

5Mi = 



-£?,M 4 



(4.19) 
(4.20) 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



(4.25) 



(4.26) 
(4.27) 
(4.28) 

(4.29) 
(4.30) 



The perturbation changes the I = 3 mass multipole, and the change is the same as for an I = 3 mass perturbation on 
a Schwarzschild background. The lower order multipoles are unchanged. 



V. PERTURBATIONS TO THE SPIN MOMENTS: SCHWARZSCHILD BACKGROUND 

When we defined mass perturbations in Sec. Mil and IIV1 we based our definition on the multipolar expansion of 
the gravitational potential in Newtonian gravity. We cannot do the same thing for spin perturbations because in 
Newtonian gravity, mass currents are not a source of the gravitational field. Instead, we create perturbations to the 
spin moments by adding an imaginary perturbation to the Ernst potential as defined in Eq. (|2.10[) : 

5 = S Schw + «Si (5.1) 

where Ss c hw = r/M — 1 is the Ernst potential for the Schwarzschild spacetime. 

When we constructed mass perturbations, we defined a set of perturbations that are proportional to spherical 
harmonic functions in the large r limit by imposing the condition that the metric satisfy the Einstein equations to 
first order. We want to define a similar set of perturbations to the spin moments. As shown by Ernst [l3|], enforcing 
the vacuum Einstein equations to first order gives the following constraint equation for Si: 

As discussed in [la ], one class of solutions to this equation includes the linearized Kerr spacetime; we discuss this in 
more detail in the appendix. We construct bumpy black holes by considering another class of solutions to Eq. (15.21) . 
Since (<9 2 Si/<9r 2 ) satisfies Laplace's equation, we can define an order I spin perturbation as one for which (c^Si/cV 2 ) 
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is an order I spherical harmonic. We can relate the perturbation to the Ernst potential to changes in the timelike 
Killing vector field by inverting Eq. (|2.10j) : 

e = §TT (5 ' 3) 

2M 2M 2 , s 

= 1 - — +*—^- ( 5 - 4 ) 

The perturbation leaves the norm of the timelike Killing vector unchanged, but it changes the twist of the timelike 
Killing vector, which was formerly zero, to 



2M 2 „ 
— r a i 



(5.5) 



The metric now has a nonzero gt^ component, which we denote o\\ 



dr \r-2MJ (r - 2M) 2 89 ' 

da 1 2M 2 s\tlO ( „ dZ 1 



( 2M\ ( 2M\ 1 

ds 2 = - ( 1 - — J dt 2 + 2cri dt d<f> + f 1 - — J dr 2 +r 2 d6 2 + r 2 sin 2 6 dcj) 2 , (5.6) 

where o\ is related to Si by 

d ( ra t \ 2M 2 sin6» 8E t 

80 - r K 2 ^-^) • (5 - 8) 
Equations (|5 . T[) and (|5.8[) overdetermine o\ ; we will use Eq. (15.71) to calculate a\ with the boundary conditions 

lim a x = , (5.9) 

r— >oo 

lim = . (5.10) 
Once we have the perturbed metric, we can calculate the multipole moments using the same procedure as in Sec. IIIII 



A. 1 — 2 spin perturbation 



We define an I = 2 spin perturbation by specifying a solution to Eq. (|5.2p that has the form of an / = 2 spherical 
harmonic: 



d 2 -! S 2 M 5 



1 



3z 2 



p 2 + z 2 



1 



dr 2 4 V 7r (p 2 + z 2 ) 3 / 2 

We can use Eq. (|3.3|) and (|3.4j) to rewrite Eq. ()5.11[) in terms of Schwarzschild coordinates: 

<9 2 Si S 2 M 15 1 



9r 2 



7rd(r,6») 3 



3(r-M) 2 cos 2 6» 

WJ) . 



Sx(r 



/J [d(r, 6>) 2 - (r - M)d(r, 0) + M 2 cos 2 9] 



M d(r, 6») 

From the perturbation to S, we can calculate the perturbation to the metric a\\ 



S 2 M 5 



, , — COS ( 
2 V 7T 



2d{r 



r-M 



d(r,6) 



r J 



(5.11) 



(5.12) 



where d(r, 8) is defined in Eq. (|3.15p . Integrating, and imposing the conditions 8a\ /dr — > and a\ — > as r — > oo 
yields 



(5.13) 



(5.14) 
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The first few multipole moments are listed below: 



SM = 0, 
SMi = 0, 



SM 2 
SM 3 

5M 5 



.1 

'4 
, 



. 



(5.15) 
(5.16) 

(5.17) 
(5.18) 

(5.19) 
(5.20) 



The I = 2 spin perturbation changes the even spin moments for I > 2, but leaves the odd spin moments and all of the 
mass moments unchanged. 



B. I = 3 spin perturbation 



Consider an / = 3 spin perturbation: 



d 2 ~! _ S 3 M 2 /7 I 
dr 2 ~ 4 



7T (p 2 + Z 2 ) 2 



5z J 



3z 



(p 2 + Z 2 f/ 2 (p 2 + Z 2 ) 1 / 2 



Then the perturbation to the Ernst potential Si is given by 



.S a 7 



(r-M)(3d(r, 9) 2 - M 2 cos 2 ' 



d(r, 6>) E 



This corresponds to a perturbation to the metric given by 



01 



S 3 M /7 



1 



2M \ r 2 (2r - 3M) - 3Mr (r - 2M) cos 2 (9 - 2M 3 cos 4 6 2r - M ~ 3M cos 2 



6 V 7T 

The first few multipole moments are listed below: 



d{r,ey 



SM = 0, 
SMi = 0, 



SM 2 = 0, 

1 

12 



SM 3 
SMi = , 



i^zS 3 M\l- 



6M 5 = 



1 

36" 



-i^S 3 M 6 \l- . 



The I = 3 spin perturbation changes only the odd Geroch-Hansen moments for I > 3. 



VI. PERTURBATIONS TO THE SPIN MOMENTS: KERR BACKGROUND 



(5.21) 



(5.22) 



(5.23) 



(5.24) 
(5.25) 
(5.26) 

(5.27) 
(5.28) 

(5.29) 



We generate a Kerr spacetime with perturbed spin moments by applying the Newman- Janis algorithm to the 
Schwarzschild spacetime with perturbed spin moments, whose metric is given in Eq. (|5.6I) . This yields 



ds 1 



2Mr iaMr > „ 
1 - ^t^cti dt 2 



E 2 



AaM r sin 



-—dr 2 + E d6 2 



i 2 i 2\2 2 a -2/1 

[r + a ) — a A sin 



- + 4(7l 

AaMr 
S - 2Mr 



(r 2 +a 2 )2Mr 



E 2 



A 



A 



tri A 



2Mr(r 



2E - 4Mr 

,2 / 



dt d<p 



sin 



d0 2 



(6.1) 
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The norm of the timelike Killing vector is 

2Mr AaMr , „. , , 

\ = l- — + ^^a 1 (r,9) . (6.2) 

As in Sec. IIV1 we cannot define the twist because the spacetime is not vacuum. For example, in the large r limit and 
expanding in a, the stress-energy tensor of the Kerr spacetime with an / = 2 spin perturbation, which we define in 
the following section, has nonzero terms 

39 4 [dcose^e 

Tee = -—aS 2 M \ = , (6.3) 

27 4 ^ cosflsin 4 ^ 
T H = -—aS 2 M \ . (6.4) 

In general, an order I spin perturbation creates nonzero terms in the stress-energy tensor that fall off like r~( l+3 \ 
As in Sec. IIV1 this allows us to define the Geroch-Hansen moments up to order Z + 1. Unlike in the case of a mass 
perturbation on a Kerr background, we will find that the twist depends on the perturbation, so we cannot go further 
in our calculation without choosing a particular perturbation. 



A. 1 = 2 spin perturbation 

We define an I = 2 spin perturbation with the perturbation potential 



S 2 M (5 
(Ji(r,6) = — y-cos< 



M 



L(r, 9, a) + M sin 2 9 L{r, 9, a) — M — d(r, 9, a) 



L(r,9,a) + M 



d(r, 9, a) 



L(r,9,a) + M 



(6.5) 



where d(r,9,a) and L(r,9,a) are defined in Eq. (|4.9I) and (|4.10[) , respectively. Now we need to construct the Ernst 
potential. The norm of the timelike Killing vector is given by Eq. (|6.2[) . As discussed in Sec. IIV1 we can write the 
curl of a nonvacuum spacetime in the form of Eq. (|4.5p . We define uj' by ignoring all terms of order in the curl. 
This yields 



UJ 



2aMcos6 S 2 M 5 1 
E ~ V 7rr 



r 2 - 2Mr + 2M 2 cos 2 9 
rd(r, 9) 



71/ 
r 



In this case, v a falls off like r 5 , so we are justified in treating uj' as the twist. The Ernst potential is 



AI 



+ t- 



S 2 M / 5 r(r - 2M) - (r - M)d(r, 9) + 2M 2 cos 2 9 



r-M-iacos6 ' " 4 V tt (r - M) 2 d(r, 9) 

Applying the procedure from the previous sections gives the multipole moments: 

SMq = 0, 
SMi = 0, 

'5 



5M 2 
SM 3 = 



i]s 2 M\ 

4 V 7T 



(6.6) 

(6.7) 

(6.8) 
(6.9) 

(6.10) 
(6.11) 



The perturbation changes the I — 2 spin moment by the same amount as an I = 2 spin perturbation on a Schwarzschild 
background, and it leaves lower order moments unchanged. 



B. I = 3 spin perturbation 



We consider an I — 3 spin perturbation: 



<ri{r,0) 



S*M 7 



L(r, 9, a) - M s 30 (r, a) + s 32 (r, a) cos 2 9 + s 34 (r, a) cos 4 



L(r, 9,a)+M 



6 



2L(r,6»,a)-M(3cos 2 6»- 1) 



d(r,9,a) 3 



L{r,9,a)+M 



3.12) 
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where 



s 30 (r,a) = 2L(r, 9, a) 3 + 3r(r - 2M) + 2M 3 
«32(r,o) = -3M(r(r - 2M) - a 2 ) , 
s 34 (r,a) - -2M 3 



3a 2 M . 



(6.13) 
(6.14) 
(6.15) 



Now we need to construct the Ernst potential. The norm of the timelike Killing vector is given by Eq. (|6.2p . As 
discussed in Sec. IIV| we can write the curl of a nonvacuum spacetime in the form of Eq. (14. 5p . We define to' by 
ignoring all terms of order aS 3 in the curl. This yields 



2aMcos9 S 3 M 2 /7cos( 



(r - M)(3r 2 - 6rM + 2M 2 cos 2 9) 
d{r, <9) 3 

In this falls off like r 6 , so we are justified in treating lu' as the twist. The Ernst potential is 



(6.16) 



M 



r — M — ia cos I 



+ i- 



S 3 M 2 17 



12 



— cos( 

7T 



3r(r - 2M)(d(r, 9) - r + M) + M 2 cos 2 9(3d(r, 9) - 2r + 2M) 



{r-M) 2 d{r,9f 



Applying the procedure from the previous sections gives the multipole moments: 

SM = 0, 
SMi = 0, 
SM 2 = 0, 

6M 3 = i±S 3 M*JI , 



SM 4 







.17) 



(6.18) 
(6.19) 
(6.20) 

(6.21) 
(6.22) 



As in the case of an I = 3 spin perturbation on a Schwarzschild background, the perturbation changes the I = 3 spin 
moment but leaves lower order moments unchanged. 



VII. CONCLUSION 

This paper extends the formal definition of bumpy black holes that was introduced in 

H and 0- We consider two 

kinds of static, axisymmetric perturbations: those that perturb the mass moments and those that perturb the spin 
moments. The mass perturbations, as defined in 0, correspond to pure multipoles in the Weyl sector. We map these 
perturbations to changes in the Geroch-Hansen moments and showed that an order I Weyl perturbation changes the 
Geroch-Hansen mass moments above Mi but leaves the lower order moments unchanged. We introduce bumps that 
perturb the spin moments, and we show that an order I spin perturbation changes the Geroch-Hansen spin moments 
above Mi. In principle, this allows us to build spacetimes whose multipoles agree with those of the Kerr spacetime 
up to some arbitrary order L but differ for I > L. 

Now that the bumpy black hole formalism has been developed, the next step is to use this framework to construct 
tests of general relativity. As discussed in [7[, changes to multipolar structure of the spacetime result in changes 
to the frequencies of geodesies. They also change the images that distant observers would see when rays are traced 
through the spacetime [l6j . Previous work has focused on the effect of perturbations to the mass moments; we expect 
that perturbations to the current moments will produce different effects. Detailed analysis will be needed to see if 
perturbed current moments also leave an observationally important imprint. For any system in which the relevant 
motions are well described by geodesies, at least on short timescales, this framework should enable us to construct 
useful and interesting tests of the black hole no-hair theorem. 
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Appendix A: Linearized Kerr 



One solution for a current-type perturbation is well known: the linearized Kerr spacetime, i.e. Kerr expanded to 
leading order in the spin parameter a. In our notation, this spacetime is given by Eq. (|5.6p with 

2aM sin 2 9 , A , 

o\ = • Al 

r 

The purpose of this appendix is to demonstrate that the framework we have developed includes the linearized Kerr 
spacetime; however, it requires a modification to the procedure we present in Sec. IVl 
To begin, note that one family of solutions to Eq. (|5.2p can be written 



—l 



Af (6) + Brf 1 (6). (A2) 



Let us choose a solution with B = and fo = cos 9. We calculate the perturbation to the spacetime u\ by enforcing 
Eq. (EU): 

do\ AAM 2 sin 9 cos 9 , A . 
69 r 



which we readily integrate to obtain 



r 

It is simple to verify that this satisfies Eq. (|5.7|) . If we now choose 



2AM 2 sin 2 9 N 
o-i = • (A4) 



A=-—, (A5) 



we identify this solution as the linearized Kerr spacetime. 
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